Confinement, Turbulence and Diffraction Catastrophes 

J.-P. Blaizot a and M.A. Nowak b 

*IPhT, CEA-Saclay, 91191 Gif-sur-Yvette, France 
M. Smoluchowski Institute of Physics, Jagiellonian University, PL-30-059 Cracow, Poland 



Abstract 

Many features of the large N c transition that occurs in the spectral density of Wilson loops 
as a function of loop area (observed recently in numerical simulations of Yang-Mills theory 
by Narayanan and Neuberger) can be captured by a simple Burgers equation used to model 
turbulence. Spectral shock waves that precede this asymptotic limit exhibit universal scaling 
with N c , with indices that can be related to Berry indices for diffraction catastrophes. 
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1. Spectral shock waves 

Recent lattice studies pQ have provided evidence for a novel universal phenomenon in 
Yang-Mills theory with a large number of colors N. For small loops, the eigenvalues of 
the unitary Wilson operator concentrate around the point z = 1 in the complex plane, 
and spread symmetrically on the unit circle towards z = — 1 as the area of the loop 
increases. At some critical area, the spectral gap closes, and the corresponding level 
density is given by a universal function (Pearcey function), whose two arguments exhibit 
scalings in N 1 / 2 and iV 3 / 4 , respectively, apparently independently on the dimensionality 
of the Yang-Mills theory. For very large loops, the level density goes towards a universal 
distribution, reflecting disorder, as suggested long ago by Durhuus and Olesen [2]. In 
a recent paper [3], we suggested a particular mechanism leading to such transition, by 
noting that in d = 2, the spectral evolution of Wilson loops is governed by a complex 
Burgers equation 

7T 

d A F + Fd e F = : F{8,A) = da p(a, A) cot((d - a) /2) (1) 

— 7T 

where the imaginary part of the resolvent F(0,A) gives spectral density p. All known 
properties of the Durhuus and Olesen transition can be easily recovered by analyzing 
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the solution of Eq. fl}, using the methods of complex characteristics, and following 
the motion of singularities in the complex plane [3]. In this description, the transition 
corresponds to the collision of two "tsunami-like" shock waves (associated to the edges 
of the spectrum) evolving symmetrically from z = 1 around the unit circle, to merge 
violently at z = — 1 at the critical area. For larger areas, the cascade of Fourier modes 
evolves in such a way that only the lowest mode, corresponding to the size of the system, 
survives, a phenomenon reminiscent of the inverse turbulent cascade. 



2. Non-linear diffusion, universality and colored catastrophes 



The complex Burgers equation describes, under certain conditions, the random walk 
of eigenvalues of large unitary matrices [4] .The role of the time t is played by the area 
of the loop. Writing the eigenvalues as z = e 1 ® , the Langevin equation corresponding to 
the random walk of unitary matrices reads 

dBi =dB l + ^-J2 cot j dt ( 2 ) 

where the first term is a Brownian noise and the second is originating from the usual 
repulsion of eigenvalues. From the Langevin equation, a Smoluchowski-Fokker-Planck 
equation is obtained using standard techniques, which, after a suitable approximation 
yields the following equation for p [5] 

d t p + de [pR[p\] = \d ee p, H[p] = -Lp.V. / dap(a, A) cot (3) 

— TV 

where P.V. denotes the principal value of the integral. Eq. is equivalent to the viscid 
Burgers equation, the last term playing the role of a viscosity term. This equation exhibits 
various scaling with N, which we now discuss. 

• Dyson gas: We rescale the density, originally normalized to N, as p = Np, and the time 
as r = Nt. Then, in the large N limit we get [5] (equivalent to Eq. [TJ 

d T p + d e [pH[p}} = 0. (4) 

Note that the viscous term has disappeared, since it was dwarfed by the effective spectral 
viscosity, v s — X/2N H"l In this regime the evolution proceeds without collisions. 

• Dyson liquid: There is however another regime, where collisions between the eigenvalues, 
taken into account by the viscous term, cannot be ignored. This happens in the vicinity 
of the shock waves, and seems to be responsible for the universal scaling at the edges of 
the spectrum. To approach this limit, we set p(6,t) = N a+1 p(x,r) with x = N a 9 and 
t = N 2a t, and get 

d T p(x, t) + Nd x (pHp) = ^d xx p. (5) 

This viscid complex Burgers equation allows us to describe the scaling properties of the 
spectral density in the vicinity of a shock. 



This value of the spectral viscosity has also been obtained by Neuberger [9] using different arguments. 
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This can be seen schematically by applying the Cole-Hopf trick to solve the equation [7] . 
The resulting spectral distribution can be written in the form 

p{x) = -^-had m ljxK(x) (6) 
iv7r 

where K{x) is the convolution of a Gaussian kernel with Kq(z) = exp{— N J z Fo(w)dw}, 
where Fq is the resolvent corresponding to the appropriate initial condition for the mi- 
croscopic spectral density p. In the vicinity of the shock waves, the resolvent can be 
approximated by a polynomial of second order (gapped phase) or third order (at the 
closure of the gap). The resulting spectral densities are therefore related to complex Airy 
or Pearcey type integrals, with well defined scaling properties. 

In fact these scaling properties appear to be universal and manifest themselves in var- 
ious physical contexts. We mention here an amusing analogy with optics. In geometric 
optics (vanishing wavelength A), the rays of light may condense on surfaces of infinite 
intensity, the so-called caustics. These correspond to stable singularities that are classi- 
fied, the lowest two corresponding to fold and cusp singularities. At finite wavelength, 
interference effects enter, and we have to use wave packets to describe the emergence of 
singularities. Such a wave packet \& is generically of the form [5] 

For the fold, (3 — 1/6 and a = 2/3, and ^> is given by an Airy function, while for 
the cusp (3 = 1/4, eri = 1/2, 02 = 3/4, and ^ is given by a Pearcey function [8]. The 
analogy with these "diffractive catastrophes" goes as follows. The limit of geometrical 
optics corresponds to the limit N — 00 (A ~ 1/^0- Rays of light correspond to the 
(complex) characteristics of the Burgers equation, whose envelopes correspond to the 
caustics and are associated to singularities of Burgers equation. The universal behavior 
in the vicinity of these singularities is captured by a Pearcey function, and the particular 
critical exponents for scaling with "angle" and "area" observed in pQ are identical to the 
Berry indices for a cusp catastrophe. 
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